Dispersion relation and spectral function of an impurity in a one-dimensional quantum liquid 
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We consider the motion of an impurity particle in a general one-dimensional quantum fluid at zero temper- 
ature. The dispersion relation fi(P) of the impurity is strongly affected by interactions with the fluid as the 
momentum approaches irnhn, iSvrfin, . . ., where n is the density. This behavior is caused by singular ±27rftn 
scattering processes and can be understood by analogy to the Kondo effect, both at strong and weak coupling, 
with the possibility of a quantum phase transition where f2'(±7rn) jumps to zero with increasing coupling. The 
low energy singularities in the impurity spectral function can be understood on the same footing. 



The study of systems in which one or a few distinguished 
degrees of freedom are simultaneously coupled to a 'bath' 
consisting of an infinite number more finds applications in all 
branches of physics. In condensed matter physics, such situa- 
tions are often called 'impurity' problems, since the canonical 
examples describe the dynamics of a defect in an otherwise 
perfect crystal lattice, or a foreign atom in a pure fluid. 

A natural distinction within the class of impurity problems 
is between the situation typical of the solid state, in which 
impurities are fixed and have infinite mass, and that found 
more often in the study of quantum liquids such as super- 
fluid ''He, where the impurity is mobile and has finite mass. 
While the former has historically provided a richer source of 
new concepts, including the related phenomena of the Kondo 
effect 11], the X-ray edge problem |i2j|3J and Anderson's or- 
thogonality catastrophe 0), the latter finds new applications 
in the physics of ultracold atomic gases. Ref. f5l provides an 
example of this new setting. A small fraction of atoms in a 
Bose condensed gas were transferred into a different hyper- 
fine state, forming a dilute system of impurities propagating 
through the condensate. With the momentum and energy of 
the impurities under independent control, a more sophisticated 
version of the same experiment could measure the spectral 
functions discussed later in this work. 

Recent experimental advances mean that the study of mo- 
bile impurities in low dimensional, strongly interacting atomic 
gases has become a possibility |6|. The purpose of this pa- 
per is to show that in one dimension such a system displays 
remarkably rich behavior The model that we will discuss ap- 
plies equally well to solid state systems such as quantum wires 
and carbon nanotubes in which the gas corresponds to a par- 
tially filled sub-band and the impurities to carriers propagating 
in an unfilled sub-band. 

The problem to be considered can be motivated by the fol- 
lowing simple calculation. Consider an impurity of mass M 
moving through a Fermi gas consisting of particles of mass ni 
and density n. With an eye to applications in ultracold physics 
the interaction between the gas and the impurity will be taken 
to be Hint = ^{'Xi ^ We take u > 0: for attractive 
interactions a bound state will form and we expect the prob- 
lem is similar but with a 'molecular' impurity. Let us find how 
the dispersion relation n^"\P) = P'^/2M of the impurity is 
modified by this interaction. At first order in u the contribu- 
tion il'^) = un is momentum independent, while the second 
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FIG. 1: Schematic dispersion relation of an impurity moving in 
K > 1 Luttinger liquid. At weak coupling a cusp persists at 
P = nn, but vanishes discontinuously at some critical coupling. 
In the strong coupling limit the dispersion takes the simple form 
S7(P) oc sin^ P/2n. ^1{P) is the ground state energy of the sys- 
tem at given P: the singular form of the impurity spectral function 
near this threshold is indicated in the inset. The dashed curve illus- 
trates the ground state in the absence of an impurity for reference. 



order contribution is (we set fi. = 1) 
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Near the Fermi points zLpp = ±7rn we have the singular be- 
havior 
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so that the derivative ^'{P) has a logarithmic singularity at 
the Fermi points. The case A/ = to is special, and we have 
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In both cases the contribution is negative. The origin of the 
singular behavior in Eq. ^ as P ^ ±Pf is the follow- 
ing. There are intermediate states in the second order cal- 
culation involving a soft particle hole pair with momentum 
~ 0, leading to a vanishing denominator. Additionally, when 
P ±Pf, a vanishing denominator arises for an interme- 
diate state in which the impurity is scattered to the opposite 
Fermi point, accompanied by a particle hole pair of momen- 
tum ~ ±2p^. We will see that this is analogous to the con- 
tribution of spin-flip processes at second order in the Kondo 
problem and corresponds to identical singular behavior as a 
function of magnetic field. 

In the above calculation, the meaning of ri(P) is clear, be- 
ing the energy of a state that goes over to the free impurity 
with momentum P at zero coupling. How should we under- 
stand ri(P) in the general case, without reference to pertur- 
bation theory? A kinematic argument shows that the ground 
state energy of the system with total momentum P may be 
identified with the impurity dispersion as long as the group 
velocity < v, the speed of sound (equal to the Fermi 

velocity w^? in this case), otherwise a lower energy state will 
be one with real particle-hole or phonon excitations. Since 
ri^°) (pp) = mvp/M the above calculation is then meaning- 
ful for M > m. By defining as the ground state energy 
at given P we obtain a periodic function il(P+27rn) = i^(P)- 

The remainder of this paper is concerned with a general 
ID quantum fluid characterized by a Luttinger parameter K, 
equal to Tm/mv for a Galilean invariant system Q. The main 
findings for ri(P) are then (see Fig.[T]): 

1. For K < 1, including the case of free fermions with 
K = 1, ri'(±7rn) = for arbitrarily weak coupling 
between the impurity and gas. 

2. For K > \ (e.g. a Bose gas), Q.'{±Tin) is finite for 
weak coupling, but jumps discontinuously to zero at 
some critical coupling. 

3. At strong coupling, when the impurity is almost opaque 
to the gas particles, the dispersion takes the simple form 
ri(P) - ri(0) cx sin^ P/2n, with the bandwidth being 
set by the amplitude for tunneling through the impurity. 

Note that the impurity propagates ballistically only at zero 
temperature, with a finite mobility diverging as cx T^^ at low 
temperatures ll8ll9llT0l. 

In addition we will find the form of the spectral function 
A+(P, w) for the addition of the impurity at momentum P. 
Since ri(P) is the ground state of the system with total mo- 
mentum P, this energy represents a threshold below which the 
spectral function vanishes. Just above the threshold the spec- 
tral function has the form of a power law with momentum- 
dependent exponent 



A+(p, uo) cx e{uo - n{p)) [lj - J7(p)]"(^^ 
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In general the calculation of the exponent a{P) is a difficult 
task. In the limit of an almost opaque impurity we obtain 



a{P) = f [f + (^) ], where is a non-universal 
parameter, equal to ttv in the Fermi or Tonks gas. Finally, 
we confirm those results that apply to the A' = 1 case for 
the Fermi gas with M — m, which is a particularly simple 
integrable system. 

We begin be developing the theory at weak coupling. As in- 
dicated above, it is convenient to describe the low energy de- 
grees of freedom of the gas as a Luttinger liquid with Hamil- 
tonian ||2l 
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with [(j){x),dy6{y)] = iTTS{x — 
nian of the system is then i/imp 



y). The total Hamilto- 

f -ffgas + up{X), where 



TJi^p = p2/2M, [P,X] -i and p{x) is the density of 
the gas. In the Luttinger liquid picture this has the form 

It is convenient to make a transformation to a frame moving 
with the impurity, which can be accomplished by the unitary 
transformation Ux — e^^^i^"^", where Pgas is the momentum 
of the gas. Then H = Hgas + (P - Pge^f /2M + up{0). 
With X now absent from the Hamiltonian P is conserved and 
corresponds to the total momentum of the system. Neglecting 
irrelevant operators arising from higher order harmonics of the 
density, the low energy Hamiltonian then takes the form 



H = 
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dx 
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In writing Eq. Q we have passed from the microscopic in- 
teraction strength u to effective couplings and ?727m, de- 
scribing forward scattering and backward scattering from the 
impurity respectively. For the (5-function interaction we have 
Ucj, — U2Tm — uat lowest order in u. We have also separated 
the zero mode contributions in (j){x) and 6{x). 
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with [00, = [^'o,-^] = i- The backward scattering term 
thus changes J by 2, corresponding to a ±27rn change in the 
momentum of the gas. If vrn — 1] < P < -rm [J + 1] 
for even integer J, the second term of Eq. (|5]l, arising 
from the kinetic energy of the impurity, favors the value 
J = J . For P = TTn{J + \), states corresponding to 
J = J tJ + 2 become degenerate in the absence of back- 
ward scattering; away from this point they have a splitting 
Ap,j^^[P-^n(:7+l)]. 

Before considering the structure of perturbation theory in 
U2Tm it is convenient to first remove the forward scattering 
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terms by a unitary transformation UlUr with ifTTI 



l^L/R = exp 



±i- 



(7) 



2ti {v =f Pj/M) 
where Pj — P — nnJ and we introduced the chiral 



modes = (p/VK ± \fKB, obeying the relations 

[0L/_R(2;),9a;'(/>^/fl(a;')] = ±27ri(5ss'(5(a; - x'). Close to a 
degeneracy point we retain only the values J = J ,J + 2, 
so that Hamiltonian is conveniently written using the pseudo- 
spin 1/2 variable — J — J — 1 so that Pj = —Trna^ and 



H = 



dx 
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where 



1 - 



and Ap J indicates that the splitting is altered from its bare 
value in a way that depends upon the short distance behavior 
of the problem and not specified by the Luttinger Hamilto- 
nian. The precise value is not important in what follows. We 

have also dropped the term (/ dx(i>dxOj in Eq. (5 i 

as it gives is an irrelevant operator that does not change our 
conclusions. 

ri(_P) for P ^ ■nn{J + 1) is just the ground state energy of 
Eq. (jsll for A p^j ^ 0. In particular, the presence of a cusp in 
ri(P j corresponds to a jump in the expectation value (cr^) oc 
ri'(P) as P passes 'Kn{J+ 1), as occurs for C/2w„ = 0. Let us 
discuss the behavior of the system as a function of U2Tm ™d 
r], before returning to our original parameters. 

A cusp is preserved for small U2T,n if the backward scatter- 
ing term is irrelevant, i.e. for > 1. To understand what 
happens for larger U2Tm, or for 77^ < 1, we observe that an 
expansion of the ground state energy of Eq. ^ U2tt7i yields 
the Anderson- Yuval expansion 1,1 2J of the Kondo problem in 
a magnetic field proportional to Ap j ||231 . For ij^ < 1 back- 
ward scattering is relevant and the ground state corresponds 
to the antiferromagnetic Kondo problem i.e. a singlet with 
((7^) = and a finite impurity susceptibility, correspond- 
ing to n\T:n{J + 1)) = 0, n"{nn{J + 1)) < 0. For 
> 1, (fj^) jumps discontinously to zero from the univer- 
sal value 77"^ as u increases past some critical value lfT3l and 
the Kondo ground state switches from ferromagnetic to anti- 
ferromagnetic. 

In terms of our original parameters, this means that for 
K < I, n'{T:n{J + 1)) = always, while for > 1 a 
transition occurs with ri'(7rn(J' + 1)) jumping to zero as u 
is increased. The determination of the critical value of u is a 
difficult problem that we do not address here. The divergence 
of Eq. (j9| as M ^ nn/v is to be expected, as in this limit 
the bare group velocity at momentum zLnn equals the sound 
velocity. This limitation does not exist for the strong coupling 
limit where the dispersion is almost flat, as we will now show. 



As u ^ 00 the impurity presents a hard wall to the parti- 
cles of the gas. If the impurity's mass were infinite, this would 
(9) correspond to the boundary conditions 0(O±) — The 
equations of motion 9 = vdx<p, 4> = vdxd then show that dxd 
vanishes as the origin is approached from either side, corre- 
sponding to vanishing current, while 9{x) is discontinuous at 
the origin: the phase difference can fluctuate wildly as the gas 
is cut in two. The value to which is pinned is determined 
by the magnitude of the forward scattering terms in Eq. (j5]l in 
the u 00 limit and is non-universal. In the simple case of a 
free Fermi gas or Tonks gas {K = 1) we have = vtt. 

In order to examine the effect of the impurity kinetic energy, 
we write the total momentum of the gas as 



-Pgas = - / dx dx9dx4> = -ne{x)\^^'l+- [ dxdxOdxC, 

J\x\>e J\x\>e 

The origin is excluded as there is a break in the fluid here. 
The second term will be unimportant at low energies, so that 
we can take the impurity term in the Hamiltonian to be 



Himp=^{P + ne{x)f,: 



The overall Hamiltonian H„ 



'imp is quadratic, 
but without discussing the solution explicitly we note 
{0{x)\q'^) Hg-^^+Hi,„p = —P/n with fluctuations about this 
value vanishing at energies < r? jM due to the pinning ef- 
fect of i/inip- Note that in this limit the ground state energy 
has no dependence on P, i.e. the impurity dispersion is flat. 

Turning now to the case of small but finite transparency, 
the term in Hamiltonian allowing tunneling of particles of 
the gas through the impurity is the Josephson-like term 
Bt = -2icos[6'(0+) - 0(0-)]. The first order correc- 
tion to the energy is then just the expectation value (iJf) ~ 

-2t(cos[6'(0+)-6'(0_)]) 2t(m/M)^/^''cosP/n, so that 

the dispersion has the form advertised earlier. 

We next discuss the behavior of the spectral function. 
The spectral function A+(P, w) is the Fourier transform of 
(0|^p(t)\E'p(0)|0), where the operator creates an impu- 
rity with momentum P, and |0) denotes the P ~ ground 
state of the gas in the absence of the impurity. Let us first 
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consider the strong coupling case. We have established that 
the impurity kinetic energy i/imp ensures that at energies 
< jM the system is described by the simple boundary con- 
ditions 



lo! = -Pl^ 

,0. KUs 



(10) 



These boundary condition can be imposed on fields that are 
continuous at the origin by the unitary transformation 



exp 



iP 



m 



TTV 
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(11) 



so that we have 

(0|vl/^(t)vl/t,(0)|0) = 0We-^"(^)*(4WZ^p(O)). 

Since {U\,{t)U p{Q)) cx \t\~^^i^) the spectral 

function has the form Eq. ^ with a{P) as given earlier. The 
result agrees with Ref. |14|, which treated the special case 
M = m,K = l (see also Ref. [15]). 

Finally, we illustrate this behavior for the case of free 
fermions coupled to an impurity with M = m. This sys- 
tem is integrable, being an extreme limit of the spin- 1/2 Fermi 
gas fT6l [T7l[T8l . As shown in Ref. |T9l, an eigenstate of a sys- 
tem of A'^ fermions can be written in the co-moving frame as a 
Slater determinant of functions fj{x) j — 1, ■ ■ ■ N satisfying 
(pj{0) — fj{L) and — <p'j{L) ~ muLpj{0) (the latter 

involving the reduced mass m/2, different from a static scat- 
terer). The (pj{x) are expanded in terms of the + 1 plane 
wave states: ipj{x) ~ X^tLo aje**^*^, where the momenta sat- 
isfy kjL — 2'iTnj — 2(5(fcj) with integers rij and phase shift 



^ . /2(fc-A) 

d(/s) = sgn(/c) + arctan 

2 V rau 



The lab frame energy and momentum are -E 2m Sj!=o 

the P — X^j^o ^J ' ^° ^'^^ spectral parameter A gives the 
former as an implicit function of the latter. In the limit L 
oo we find £;(A) p\L/{ioTTm) + f^(A) with 
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We find the following i) Eq. (j2|i holds for 



Pf _ iPp 
2m ZTrm' 



- cos2 {t:P/2pf) for 



0, ii) Vl{P) = 
0, and 
Thus 



oo ox pp —> 

iii) ^(P) - ir - as P ^ ±pp 

^l'{ztpp) — for arbitrary u. The fact that il{±pp) — 
Pp/2m for all it is a consequence of SU{2) symmetry. 

The extremely simple form of the eigenstates allows us to 
extract the behavior of the spectral function easily. Indeed 
it is possible to choose the Lpj {x) so that they approach plane 
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FIG. 2: Dispersion relation (top) and spectral function exponent 
(bottom) for the exactly soluble limit of an impurity of mass 
M = m moving in a Fermi gas, corresponding to Luttinger 
parameter K = 1. Results are given for 7 = mu/h^n = 
0.01(0), 0.2(n), 0.5(A), 1(0), and 2(*). The dashed curve gives 
the exponent a{P) — 



[wk) ] ^^li'l ill 7 — > 00 limit. 



wave states with forward scattering phases 6{k) in the L ^ 00 
limit |fT9l . The form of the spectral function exponent there- 
fore coincides with the classic result for the deep hole cor- 
relation function in the X-ray edge problem 13]: a{P) = 
[S{pp)^ + S{—pp)^) /vr^. In the u ~^ 00 limit one finds 

Si±pp) = 1^ T f and thus a(P) = i[l + ( Ji)'] in agree- 
ment with our general result. The form of il{P) and a{P) 
are shown in Fig. |2] for various coupling strengths. The per- 
sistence of edge singularities for a massive impurity in one 
dimension was noted some time ago for P = Ii20i i21il . 

In summary, we have shown that the motion of an impu- 
rity in a one dimensional quantum gas has some interesting 
features that have their origin in singular backward scatter- 
ing processes at particular momenta. This system has been 
recently realized in an ultracold atomic gas ll22l . 
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bekov, Alex Kamenev, and Leonid Glazman. 
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